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Abstract: We present a Monte-Carlo approach to prompt-photon production, where photons 
II ■ and QCD partons are treated democratically. The photon fragmentation function 

Q-ij is modelled by an interleaved QCD+QED parton shower. This known technique 

is improved by including higher-order real-emission matrix elements. To this end, 
we extend a recently proposed algorithm for merging matrix elements and trun- 



cated parton showers. We exemplify the quality of the Monte-Carlo predictions 
by comparing them to measurements of the photon fragmentation function at 
LEP and to measurements of prompt photon and diphoton production from the 
| Tevatron experiments. 
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O ■ 1 Introduction 



The measurement of final states containing photons at large transverse momenta plays a key role in collider 
experiments. Most prominently at hadron colliders inclusive diphoton or diphoton + jet signatures are 
promising channels to search for a light Higgs boson [T] . Signatures with photons might furthermore provide 
access to physics beyond the Standard Model like supersymmetry or extra spatial dimensions [TJ. Less 
spectacular but extremely important though, photon + jet final state can be used to determine the absolute 
energy scale of \ow-px jets and to constrain the gluon distribution inside the beam hadron |14j . The 
success of the outlined physics menue however strongly depends on our ability to thoroughly understand and 
accurately simulate such prompt-photon production processes in the context of the Standard Model. 

In the framework of perturbation theory the mechanism of hard-photon production is twofold. A photon can 
be well-separated from any other particle in the collision, which makes it possible to describe the reaction 
with fixed-order matrix elements. The fact that these matrix elements include initial- and/or final-state 
QCD partons necessitates an all-orders resummation of large logarithmic QCD corrections, which are then 
absorbed into Parton Distribution Functions (PDFs) and Fragmentation Functions (FFs). Due to its van- 
ishing mass, a photon can also be infinitely close to an initial- or final-state QCD parton. The related 
singularities in hard matrix elements are absorbed into process-independent photon fragmentation func- 
tions [17] , describing the transition of a QCD parton into a bunch of hadrons and a not well-separated 
photon during the process of hadronisation. Due to the nonperturbative nature of the hadronisation pro- 



cess, parton-to-photon fragmentation functions contain a nonperturbative component and must therefore 
be determined from experimental data. Their evolution with the factorisation scale hf,~/ can however be 
calculated perturbatively. While the description of hard photons through matrix elements is said to yield 
the direct component of photon observables, the description by fragmentation functions gives the so-called 
fragmentation component. Both components are related by factorisation and must be combined to obtain a 
meaningful prediction of QCD-associated photon production. 

The standard method to theoretically devise a meaningful prompt-photon cross section is to reflect certain 
experimental photon-isolation criteria in perturbative calculations. However, one must allow for a minimal 
hadronic activity in the vicinity of the photon. Only then it can then be ensured that all QCD infrared 
divergences are properly cancelled. Several such criteria are on disposal, e.g. the cone approach [T51 [T3] . 
the democratic approach [50] and the smooth isolation procedure [33]. For both, single- and diphoton 
production at hadron colliders, the complete next-to-leading order QCD corrections to respective direct and 
fragmentation components are known |18025ff2"8] . The parton- level Monte Carlo programs JETPHOX [35] and 
DlPHOX [32] implement these NLO results numerically and allow the user to choose from different photon- 
isolation criteria. The NLO corrections to the direct channel 77 + ljet have been calculated in The 
results for the loop-induced processes gg — ¥ 77 [35j and gg — > 77c/ [36j are also available. Beyond calculations 
at fixed order in the strong coupling large efforts were spent on the evaluation of soft-gluon emission effects 
and on the resummation of corresponding large logarithms. Soft-gluon resummation up to next-to-next-to 
leading logarithmic accuracy is taken into account in the program RESBOS [39]. The analytic result for 
resumming threshold logarithms was presented in |42) while small- a; logarithms have been studied in |44) . 
Only recently a first study of the prompt-photon process in the framework of Soft-Collinear Effective Theory 
has been presented [35] . 

Let us note, that there is a further source of final-state photons, namely decays of hadrons, such as it or 77. 
However, such non-prompt production processes can to some approximation be separated from the other two 
experimentally and measurements are usually corrected for these effects. This is the case for all experimental 
data referenced in this work. 

In this publication we pursue a different strategy of simulating final states including photons. We account for 
the hard-production process, the QCD evolution of initial- and final-state partons, as well as the transition 
of QCD partons into hadrons by means of a multi-purpose Monte-Carlo event generator. In this context the 
lowest-order matrix elements for single- and diphoton production supplemented with QCD parton showers 
correspond to the above mentioned direct component. The fragmentation contribution is modelled by the 
incorporation of QED effects into the parton shower. In fact, a generic algorithm to treat photon radiation is 
also given by the approach of Yennie, Frautschi and Suura [46] ■ This scheme is particularly suited to compute 
logarithmic corrections arising from soft photon radiation, where the coherent emission off all QED charges 
involved in the process plays an important role. In this publication, however, we are primarily interested 
in the production of hard, well-separated photons. Such emissions need to be treated by an improved 
algorithm, see for example [37] . We therefore choose to simulate photon radiation using a dipole-like QED 
shower model. This approach only presents a primitive approximation to soft photon effects, but is easily 
realised and no additional free parameters arc introduced in the parton-shower algorithm, cf. |48j . Similar 
methods are employed in most contemporary shower programs |50j . An apparent advantage is, that this 
method also allows for a direct comparison with experimental data since it yields predictions at the level of the 
experimentally observed particles. In particular the parton-to-photon fragmentation functions are explicitly 
modelled this way. As a consequence no further corrections accounting for the non-perturbative parton-to- 
hadron transition need to be applied and again no additional free parameters need to be introduced. This 
is crucial also for the validation of a separation of non-prompt photons from prompt photons as mentioned 
above. The democratic treatment of partons and photons in this approach combines the direct and the 
fragmentation component in a very natural way. It is well suited for comparison to experiments, where it is 
often necessary to study the impact of photon-isolation criteria. 

An apparent disadvantage of the approach is that it relies on lowest-order matrix elements only and corre- 
spondingly higher-order QCD corrections are taken into account in the approximation of the parton shower 
only. We improve on this deficiency by including higher-order real-emission matrix elements. Parton-shower 
simulations supplemented with multi-leg matrix elements have become a standard tool for the description of 
QCD radiation accompanying the production of massive gauge bosons [S3] or coloured heavy states [SpJ. In 
this line we extend the formalism originally presented in |62] to the case of prompt-photon production. This 
process, as stated before, introduces the additional complication of a second source of photon production - 
the fragmentation component - which is not present for massive gauge bosons. While in most cases of W- 
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or Z-boson production the massive boson is the hardest object in the interaction, the photon is unlikely to 
play this role in most prompt-photon events. We lay out a formalism that is capable of coping with this 
situation and allows to consistently combine tree-level matrix elements of variable photon and QCD parton 
multiplicity with a combined QCD+QED parton-shower model. Photons and QCD partons are treated fully 
democratically in this scheme. 

The outline of this paper is as follows. In Section[5]we introduce our method for simulating photon production 
using a dipole-like parton shower model. Section [3] presents the formalism for combining QCD+QED matrix- 
elements of different final-state multiplicity with the parton shower. We also discuss how to efficiently 
incorporate a given photon-isolation criterion. In Section [4] we present the results of our Monte Carlo 
analysis and discuss the interplay between the direct photon contribution and the fragmentation component 
in our approach. Finally, Section [5] contains our conclusions. 



2 Interleaved QCD+QED parton evolution 

In this section we briefly recall a formalism which is used in our simulations to generate the combined 
QCD+QED parton evolution. Similar approaches are implemented in most contemporary shower pro- 
grams [50] . For simplicity of the argument we focus on pure final-state evolution. Note that any parton 
shower algorithm is uniquely defined by three ingredients: The first is the Sudakov form factor, A a (/i 2 , Q 2 ), 
i.e. the probability for a given parton, a, not to radiate another parton between the two evolution scales 
Q 2 and . The second is the ordering or evolution variable. While the Sudakov form factor defines the 
anomalous-dimension matrix, and thereby the functional form of logarithms which are resummed, the evolu- 
tion variable selects their argument, i.e. it defines a "direction" in the phase space, into which the evolution 
is performed. The third ingredient of a parton-shower model is the method, which is applied in order to 
reshuffle momenta of already existing partons when one of them goes off-shell to allow for a branching 
process. 

For QCD parton evolution we choose to employ the parton-shower algorithm of [63], including the ordering 
variables for initial-state evolution proposed in [6 2) . This means that the Sudakov form factor for final-state 
evolution reads 

[ J ^o -L Jz - i,k 

where 

K?™. (5, ki) = j( k 2 ~ z ) ( V Q V D (k 2 2) ) and ~ z = P*P\ . (2 ) 

(i 3 )i,k\ > -u 2jT -L' " (Pi+Pj)pk 

The Jacobian factor J and the spin-averaged dipole functions (V) are defined in 63J . The sums run over all 
possible splitting products i and all possible spectator partons k of the splitting parton (ij). The ordering 
parameter is the invariant transverse momentum squared 

ki = (Q 2 -m 2 ~m 2 - m 2 k ) y ij)k - Zi) - (1 - z^ 2 m 2 - z 2 m 2 , (3) 

where Q = pt + pj + pk, y%- h k = PiPj/ {PiPj + PjPk + PkPi) and m are the parton masses. 

Since QCD and QED emissions do not interfere, their corresponding emission probabilities factorise trivially. 

A combined QCD+QED evolution scheme is thus obtained by employing the combined Sudakov form factor 

A(vl Q 2 ) = AQ CD ( M 2 , Q 2 ) AQ ED ( M 2 , Q 2 ) , (4) 

where 

Kg^ki) = ^ J(ki,i)(V ( ^(ki,z)) . (5) 

Note that we use spin-averaged dipole functions, not only in the QCD, but also in the QED case. One 
possible improvement of the present algorithm would therefore be to include the spin-dependent splitting 
kernels. However, in the domain of hard-photon radiation that we are interested in, this can simply be done 
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by employing the full real-emission matrix element instead. No special treatment is therefore necessary in 
the parton shower. 

The functional form of the spin-averaged splitting kernels is largely constrained by the infrared singularity 
structure of one-loop QCD amplitudes. It is, however, not fixed and one has the freedom to incorporate non- 
singular pieces, which can help to improve the predictions of dipole-shower simulations, cf. |64j . Likewise, 
the construction of the splitting kinematics is largely constrained by the phase-space variables selected in 
the splitting. It is, however, not fixed and one has one additional degree of freedom, which corresponds 
to the choice of the angular orientation of the splitter-spectator system with respect to the remaining 
particles. One of the most prominent criticisms of the new dipole-like parton-shower models is the seemingly 
unphysical recoil strategy employed in configurations with initial-state splitter and final-state spectator. 
This strategy is entirely due to a choice for the momentum mapping between leading-order and real-emission 
kinematics, which was initially proposed in [65] . The transverse momentum of the emitted parton is thereby 
compensated by the spectator, leaving not only the virtuality of the splitter-spectator system invariant, but 
also its complete four-momentum. One can easily imagine a different recoil strategy, where the transverse 
momentum of the emission is instead compensated by the set of all final state particles. Such an approach 
was recently suggested in [57] for the case of massless partons. We extend it to the fully massive case in 
Appendix lAl and investigate the corresponding effects on the Monte-Carlo results in Sec. HI 

It is otherwise straightforward to extend the above algorithm to initial-state showering. The only subtlety 
in this context arises from the fact that the fully democratic approach pursued here also allows initial-state 
photon splitting into a quark-antiquark pair, with the quark (antiquark) initiating the hard scattering. In 
this case parton distributions which incorporate QED effects are in principle necessitated. Even though such 
PDF fits exist (e.g. [55]), the corresponding effects on physical observables should be very small, such that 
the usage of PDF's without QED contribution does not pose a conceptual problem. 

An apparent disadvantage of the above algorithm for generating QED emissions using a parton shower is 
the low efficiency with which isolated photons will be produced. This problem is dealt with in Appendix iBl 
where we introduce a method to enhance the corresponding emission probability, at the price of generating 
weighted events. 



3 Merging QCD+QED matrix elements and truncated showers 



In this section we discuss an algorithm for merging tree-level QCD+QED matrix elements and parton 
showers, based on the method proposed in (62]. We treat photons and QCD partons democratically, i.e. 
higher-order tree-level matrix elements can be of order a n a™ compared to the leading order. If n > 0, they 
may contribute to an observed hard-photon final state. In this respect, the inclusion of higher-order real 
corrections corresponds to shifting the simulation of hard-photon production from the parton-shower to the 
matrix-element domain. 

The merging approach presented in |62] is essentially based on replacing the splitting kernels of the parton 
shower by the appropriate ratio of full tree-level matrix elements in the domain of hard-parton radiation. 
This domain is identified by simple phase-space slicing. The slicing parameter, the so-called jet criterion, is 
given in terms of parton momenta Pi, pj and pk as 

/<„'/, 

2 



2piPj min 



kJQ c k -+C k - 



where Cf 



(Pi 



Pk)Pj 
1 



2p l p J 



if 3 = 9 



(6) 



else 



For initial-state partons, one considers the splitting process a — > 
momentum of the combined particle (aj) given by p a j — p a — Pj, 



(aj) j instead of (ij) 
one then defines 



ij. With the 
C, fe s - . The 

minimum in Eq. © is over all possible colour partners k of the combined parton (ij). The jet criterion Q?- 
identifies - to leading-logarithmic accuracy - the most likely splitting in a dipole-like parton cascade leading 
to the set of final-state momenta {p}. The phase-space slicing is now implemented by selecting a cut value 
Qcut and defining the evolution kernels JC ME and JC PS for matrix-element and parton-shower domain as 



^fiiKfc(^ k l) = fc(ij)i,k(z,ki) 6 Q cut - Qy(5,kj_) 



cu I 
2 



(7) 
(8) 
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The kernel JC is then replaced by appropriate ratios of tree-level matrix elements up to a given maximum 
multiplicity. A detailed description of the corresponding algorithm can be found in 62 . It is obvious that 
the same procedure can be applied to QED emissions, once they are resummed by the parton shower using 
Eq. dU). It is then in principle possible to define two separate slicing cuts, and Q^ D , which account 

for the merging of QCD and QED tree-level matrix elements with the parton shower, respectively. Within 
the context of this work, we choose to leave these slicing cuts identical, since the typical "hardness" of a 
hard well-separated final-state photon is similar to the one of a final-state QCD jet. 

In prompt photon production processes we might be confronted with a situation which cannot arise in pure 
QCD events, namely that a single, perturbatively produced particle - the photon - is identified out of 
potentially many particles forming a broad jet. Several methods exist to achieve this identification. In the 
democratic approach |20) final state particles are clustered into jets, treating photons and hadrons equally. 
The obtained object is called a photon or a photon jet, if the energy fraction z — E 1 /(E 1 + Ehad) of an 
observed photon inside the jet is larger than an experimentally defined value z cu t. In the cone approach [T5in"5] 
photons are required to have a minimum transverse momentum and to be isolated from any significant 
hadronic activity within a cone in ry-0-space. Minimal hadronic activity in the vicinity of the photon (adding 
of the order of a few GeV to the total transverse momentum in the cone) must thereby be admitted to ensure 
the infrared finiteness of observables. 

While the jet criterion Eq. ([6]) works very well also for photons defined by the democratic approach, in 
the case of the cone approach it might not be appropriate to separate matrix-element and parton-shower 
domain. Note that the main idea of the merging procedure is to improve the parton-shower prediction with 
fixed-order matrix elements in those regions of phase space which are relevant for the analysis of multi- 
jet (multi-photon) topologies. In this respect, it is certainly desirable that experimental requirements are 
reflected by Qf j. This is possible because the jet criterion. Eq. (ii). is not fixed, but rather chosen conveniently 
to reflect the singularity structure of next-to-leading order real-emission amplitudes in QCD [62] . Moreover 
it is a flavour- dependent measure, which allows to redefine it just for branching processes involving photons. 
The most common experimental requirements of a minimum transverse momentum and an isolation cone in 
77-0-space could for example be reflected by 

Q% = mm^l-^} hl D2 ^ and Ql = Pl,i , (9) 

where the first equation applies to final-state photons and charged final state particles, while the second 
applies to photons and charged beams. Note that Eq. © is essentially equivalent to a longitudinally invariant 
jet measure |69) . One can now increase the ratio of photons produced through matrix elements over photons 
produced in the shower by simply lowering the value of Q^ D . A convenient way to obtain the largest 
fraction of events from hard matrix elements is to require a jet separation below the experimental cut on 
the photon transverse momentum and by setting D lower than the radius of the experimentally imposed 
isolation cone. 

4 Results 

In this section we apply the event-generation techniques introduced above to prompt-photon production 
at lepton and hadron colliders. In Sec. 14.11 we study the capability of our proposed QCD+QED shower 
algorithm to reproduce the scale dependent photon fragmentation function measured in hadronic Z° decays 
at LEP. In Sec. l4.2l and Sec. 14.31 we turn the discussion on single- and diphoton final states at hadron colliders, 
respectively. Besides quantifying the size of the different core-process components in the democratic shower 
approach we elaborate on the impact of real-emission matrix-element corrections incorporated using the 
formalism described in Sec. [3] All results presented here are obtained using the Monte-Carlo event generator 
SHERPA [7T] in a setup described in Appendix [C] 

4.1 The photon fragmentation function 

A crucial benchmark for the combined QCD+QED shower algorithm introduced in Sec. [5] is posed by the 
requirement to reproduce the scale-dependent photon fragmentation function -D 7 (z 7 ,y C ut) [10], where z 1 is 
the photon's energy fraction with respect to its containing jet and y cu t a resolution scale, given e.g. in the 
Durham scheme. This observable was measured to very high precision in hadronic Z° decays by the ALEPH 
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collaboration [73] ■ In this analysis events are selected where all final-state particles are democratically taken 
into account for jet finding. The events are subdivided into 2-jet, 3-jet and > 4-jet topologies with at least 
one reconstructed jet containing a photon where the photon carries at least 70% of the jet energy (z 7 > 0.7) 
and Ery > 5 GeV. The resolution measure y cut is varied between 0.01 and 0.33. The measured data is 
statistically corrected for residual hadronic decay backgrounds and initial-state radiation off the incoming 
lcptons. 

Figure [1] shows a comparison between our hadron- level Monte-Carlo results and the data from [75]. In 
the left column the z 7 distribution for 2-jet events at four different y cu t values, namely 0.01, 0.06, 0.1 and 
0.33, are shown. The right column shows corresponding results for 3-jet events at y cu t — 0.01, 0.06, 0.1 
and > 4-jet events at y cut = 0.01. For all the data z 1 — \ corresponds to completely isolated photons, 
which is reflected by a strong peak in the z 7 distribution. At the parton level the lowest-order contribution 
to fully isolated photon production corresponds to qqj final states where the quarks form one jet and the 
photon makes up the other one. At the hadron level however this sharp peak gets somewhat broadened by 
hadronisation effects due to the association of soft hadrons with the photon by the jet-clustering algorithm. 
Our Monte-Carlo simulation agrees very well with the data for the measured z 7 range in the three topology 
classes at the given jet resolutions. This can be seen as a strong indication that the proposed QCD+QED 
shower scheme is indeed appropriate to describe hard photon radiation. 

The above analysis is especially tailored to study the fragmentation component of the prompt photon pro- 
duction mechanism. The key point in this respect is the application of a democratic jet-clustering procedure 
which is in one-to-one correspondence with the democratic approach used to compute the photon-production 
rate from theory [5D]. It is therefore not obvious, that the democratic approach, also underlying our simu- 
lation, performs well in experimental situations where the photon must pass a strict isolation criterion, as 
discussed in Sec. [3] We will now turn to investigate such isolated photon final-states in more detail, focusing 
on their emergence at hadron colliders. 

4.2 Prompt-photon hadroproduction 

The inclusive production of isolated photons has been measured over a wide range of photon transverse 
energies by the CDF and D0 experiments at the Fermilab Tevatron at yfs = 1.96 GeV. In [73] CDF has 
presented a measurement covering |?7 7 | < 1.0 and transverse energies between 30 < £t,7 < 400 GeV. The 
photon isolation criterion used corresponds to the requirement that the additional transverse energy found 
in a cone of R = -J (A0) 2 + (Arj) 2 = 0.4 around the photon is less than 2 GeV. A similar D0 measurement 
was described in [75]. It covers photons of transverse momentum px n > 23 GeV up to pr.-y = 300 GeV 
and \rjj\ < 0.9. Photon isolation is implemented by demanding (-Er = o.4 — -^=0.2)/ Er—o.2 < 0.1, where Er 
is the total energy found in a cone of size R around the photon. Both measurements have been corrected 
to particle level and the dominant background of photon production from hadron decays, such as tt° 77 
and 77 — > 77, has been subtracted. We therefore attempt a comparison with Monte-Carlo predictions at the 
parton level after jet evolution. 

Figure [H compares the data for d 2 a/(dET,- y dr] 1 ) from [71], respectively d 2 a/(dpT,jdri 1 ) from [75], to our 
parton-level Monte-Carlo results, obtained using leading-order matrix elements in the democratic approach 
combined with QCD+QED shower evolution. In addition to the total result (red histograms), we display 
contributions from the different classes of partonic core processes, i.e. from dijet production (jj — > jj), single- 
photon production (jj jj) and diphoton production (jj — > 77). Taking into account the uncertainties 
of the measurements and the finite Monte-Carlo statistics in the high-i£r~ bins our calculation agrees well 
with the data. For the CDF measurement the data has a somewhat steeper slope at low £t )7 and the 
Monte-Carlo calculation pronounces the high £t,7 end of the spectrum. Regarding the different sources 
of final-state photons in our theoretical model, the main contribution to this observable stems from single 
photon production. But even though strict isolation criteria are applied, there is a considerable fraction 
of dijet events, where a hard, isolated photon is produced during the parton-shower evolution in both data 
samples. This substantiates the argument that the combined shower scheme is crucial for a proper description 
of such photon final states. The diphoton core process on the other hand is negligible here. 
We now turn to study the impact of higher-order real-emission matrix elements on the results. Therefore we 
supplement the pure parton-shower evolution by tree-level matrix elements with up to two additional light 
QCD partons or photons using the matrix-element parton-shower merging formalism described in Sec. [3] 
The comparison with measurements from CDF and D0 is shown in Fig. [3] Besides the total results (red 
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histograms) we again present the sub-contributions assigned to matrix-element core processes with exclusively 
0, 1 and 2 photons plus a variable number of additional QCD partons. When comparing to Fig. [3J where the 
pure shower result is shown, it is apparent that the majority of events with a dijet core process in the shower 
simulation is now ascribed to matrix-element cores with one or two photons plus additional QCD partons. 
Thereby, what might traditionally be called fragmentation component gets significantly reduced and turned 
into what is called direct component. For this very inclusive measurement we observe no strong variation 
of the total result due to the inclusion of real-emission matrix elements. The biggest effect is a somewhat 
larger inclusive rate for the merged samples < 20%. The shape of the distributions is preserved. This in fact 
has to be understood as a highly non-trivial consistency check of our merging formalism. 

4.3 Prompt-diphoton production 

An interesting further testbed for the democratic merging approach is diphoton production at hadron col- 
liders. The CDF collaboration has measured properties of the corresponding final states in some detail. The 
analysis presented in |76] selects leading/subleading photons with transverse momenta larger than 14/13 
GeV. Those must be isolated from any significant hadronic activity within a distance of R = 0.4, by requir- 
ing the hadronic activity within this cone to yield Et < 1 GeV. For the selected events the invariant mass 
and transverse momentum of the photon pair are analysed as well as the azimuthal separation between the 
photons. 

It is worth noting that our MC simulations include the loop-induced contribution gg — > 77. It has been 
shown |76] that its main influence is seen in the invariant mass spectrum around 30 GeV where it accounts 
for a significant enhancement of the cross section. 

To again exemplify the importance of the fragmentation contribution even for the required isolated photons in 
this analysis, Figure E] compares our Monte-Carlo prediction for leading-order matrix elements plus shower 
evolution. It displays the contributions from the previously introduced classes of matrix elements (i.e. 
33 33> 33 13 an d jj — > 77). It is evident, that the democratic treatment of photons and QCD partons 
is absolutely mandatory to describe these observables. 

From a theoretical perspective this reaction is interesting because the diphoton system does not have a 
transverse momentum when described by leading-order matrix elements. Hence its transverse-momentum 
spectrum depends strongly on the proper inclusion of higher-order effects. In addition, the azimuthal angle 
gives a measure for the correlation of the two photons which is also sensitive to higher-order corrections. 
Especially in the region of large transverse momenta or large decorrelation one expects these corrections to 
be better described by matrix elements than by the parton shower. 

In this context the parton-shower kinematics might become important, because the recoil scheme discussed in 
Sec. Splays an important role for the generation of transverse momentum for the diphoton system. Thus, as 
a first step, Figure [5a| compares parton-level Monte-Carlo predictions using two different splitting kinematics. 
We observe that both, the algorithm outlined in Appendix[21 denoted "Scheme 1" , and the method proposed 
in [53] , denoted "Scheme 2" , have difficulties describing the critical regions mentioned above. 
We show in Fig. I5bl that with the inclusion of higher-order real-emission matrix elements, the simulation 
is able to describe the measurement significantly better. Especially the transverse-momentum distribution 
exemplifies two unique features: The resummation of large logarithms correctly reproduces the Sudakov- 
shape of the low-j?j_ region which is not possible with fixed-order calculations. At the same time exact matrix 
elements allow for a consistent simulation of the high-px tail where a traditional parton-shower approach 
would fail. Also the decorrelation between the photons is now matched very welfl. The simulation becomes 
largely independent of the precise implementation of the parton-shower kinematics, an effect which is also 
observed in |64) . This happens because in the relevant part of the phase space, hard matrix elements are 
employed to define the kinematics of the diphoton system. In this way the merging algorithm of Sec. [3] 
can be used to eliminate theoretical uncertainties in the parton-shower model employed. At the same time 
we substantiate our introductory statement, that within the domain of hard-photon production the parton 
shower can easily be corrected using higher-order matrix elements. 

lr The deviations in the second and third bin might be due to large statistical fluctuations in the measurement and seem to 
be contradicted by an earlier albeit unfortunately unpublished D0 measurement I77| . 
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Figure 1: The z 7 distribution measured in hadronic Z° decays by ALEPH [73 for 2-jet, 3-jet and > 4-jet 
events at different Durham resolution y cut . The theory result corresponds to QCD+QED shower 
evolution of the leading-order qq process, taking into account hadronisation corrections. 
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Figure 2: Inclusive photon transverse energy distribution compared to data from CDF [73] (left) and 
D0 [75] (right). The contributions of different classes of leading-order core processes are also 
displayed. For the notation used cf. the main text. 
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Figure 3: The inclusive photon transverse energy obtained with a QCD+QED shower simulation supple- 
mented by real-emission matrix elements with up to two additional QCD partons or photons 
(denoted 2 -> 2,3,4) is compared to data from CDF [73] (left) and D0 [75] (right). 
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Figure 4: Properties of diphoton events measured by CDF [7B]. Displayed are the sub-contributions from 
different leading-order matrix elements and their sum. 



5 Conclusions 

We have presented Monte-Carlo algorithms for the precise simulation of hard-photon production in collider 
experiments. Using interleaved QCD+QED evolution in a dipole-like parton shower enables us to simulate 
the photon fragmentation function within a general-purpose Monte-Carlo event generator. Comparison with 
data from the ALEPH experiment exemplifies the quality of the approach. 

To diminish intrinsic uncertainties of parton-shower models, we have supplemented our simulation with 
higher-order tree-level matrix elements. To do so, an existing algorithm for QCD has been extended to 
include QED democratically This can be seen as a first step towards a unified prescription for treating 
strong and electroweak radiative corrections. At the same time it provides a natural way to simulate 
hard-photon production, where the fragmentation component is described consistently by the combined 
QCD+QED resummation in the parton shower supplemented with a non-perturbative hadronisation model. 
We have employed this procedure to analyse prompt-photon production in hadron-hadron collisions and find 
an improved description of Tevatron measurements. Due to a much larger phase space available for radiative 
corrections, these effects should become even more pronounced at LHC energies. 



Acknowledgements 

We would like to thank Thomas Binoth, Gudrun Heinrich, Thomas Gehrmann and Frank Krauss for numer- 
ous fruitful discussions. Help from Thomas Binoth with the code for gg — > 77 is gratefully acknowledged. 
We are grateful to Frank Krauss and Thomas Gehrmann for valuable comments on the manuscript. We also 
thank Dmitry Bandurin and Michael Begel for their clarifications about the D0 measurement. FS would 
like to thank the ITP Zurich for its hospitality while this work was finalised. SH acknowledges funding 
by the Swiss National Science Foundation (SNF, contract number 200020-126691) and by the University of 
Zurich (Forschungskredit number 57183003). The work of FS was supported by the MCnet Marie Curie 
Research Training Network (contract number MRTN-CT-2006-035606). SS acknowledges financial support 
by BMBF. 



A Dipole-splitting kinematics 

In this appendix we derive alternative splitting kinematics for dipole-like parton showers in the spirit of |63j . 
In addition to the massless case proposed in [6 7) we give explicit formulae for the fully massive case, which 
plays an important role for truncated-shower algorithms and in processes involving heavy quarks |78j . Ef- 
fectively, only one dipole configuration is considered, i.e. branching final-state partons with the spectator 
parton being in the final state. In this case we closely follow the proposal in [63] . which is inspired by the 
original kinematics of the massive Catani-Seymour dipoles [SS]. Kinematic relations for all other dipole 
configurations are then derived using crossing relations. 
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Figure 5: Properties of diphoton events measured by the CDF collaboration [76]. Figure (a) compares the 



influence of different parton-shower kinematics when using leading-order matrix elements. Fig- 
ure |(b)| shows the same comparison for merged event samples with up to two additional particles 
in the matrix element-final state. Scheme 1 refers to the algorithm outlined in Appendix ^ 
scheme 2 stands for the original implementation |63j . 
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Q Q 
Figure 6: Schematic view of the splitting of a final-state parton with a final-state specta- 
tor. The blob denotes the hard process from which the three partons emerge. 
Since the four-momentum of the splitter-spectator system is conserved, the 
pair can be regarded as decay products in a two-particle transition of a virtual 
particle with momentum Q = p^ + pk ■ 



We consider the process depicted in Fig. [SI where a parton ij, accompanied by a spectator parton k, splits into 
partons i and j, with the recoil absorbed by the spectator k. We define the combined momenta p.y = pi +pj 
and Q — pij + pk and the variables 



Vij-k 



PiPj 



PijPk 

Thus, we immediately obtain 

_ Vij,k 



PiPk 
PijPk 



1 + Vij,k 



(Q 



(mf+m 2 ) 



Now let the light-like helper vectors I and n be given by (cf. [75] ) 



Pij - OlijPk 



Pk - a k Pij 



1 



aijCtk 



1 



OtijOtk 



where 



and 



lij,k 



= 2 In = 



(Q 2 



I) + sgn (Q 2 - Sij - m 2 k ) ^\(Q 2 , SlJ ,m 2 k ) 



with A denoting the Kallen function A(a, b, c) = (a — b — c) 2 — 4 be. 

The momenta pi and pj can then be expressed in terms of I, n and a transverse component k± as 



ot =Zil» + 



2Vn 



(1 - Zi) I" + 



i 2 + k 2 x 



2lr 



k^ , 



(10) 



(11) 



(12) 



(13) 



(14) 



The parameters z% and of this decomposition are given by 



Q 2 



mi 



^\(Q 2 ,Sij,m 2 k ) 
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ki = (Q 2 - 



2 2 2\ 



\Hj,k 
Vij^k 
+ Vij.k 



Vij,k + 



2mf 



Q 2 



V ""k 
2^2 



Zi (1 - Zi) - (1 - zt) m 



2 2 

z t rrij 



(15) 



Equations (fTTj) and (| 15[) are valid for all dipole configurations, i.e. initial and final-state branchings with the 
recoil partner being either in the initial or in the final state. The corresponding mapping of variables t]ij k 
and ^i,jk onto those defined for massless partons in [SS] is listed in Tab. [TJ Note, that in the massless case 
these kinematic relations reproduce the results of (67) . 
Constructing a shower emission proceeds as follows 



1. Determine T)ij t k and £i,jfc from evolution and splitting variable. 

Compute Sij according to Eq. (TTTj) and z% and k^_ according to Eqs. (fT5]l . 
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Configuration 




Vij.k 


FF 


Zi 


Vij,k 

1 - Vij,k 


FI 


- Zi 





Configuration 



IF 



II 



£,j,ak 



Vja,k 



1 _ 1-Uj 



%jk,a u>j 



%j k , a 



i 



%j,ab i Vj 



Table 1: Mapping of variables for Eqs. (HTJ and (lT5j) . 



2. Construct pj, according tc0 [65 



pk - r ^ — Q va(o>r) + — 2^ — ■ (16) 



3. From and construct the light-like momenta I and n according to Eqs. (1121) . 
Determine the new momenta pi and Pj according to Eqs. f) 1 4[) . 

Due to the possibly vanishing denominator of £ijk an d rjij^ in the case of initial-state splitters with final- 
state spectator, the corresponding kinematic relations occasionally become numerically unstable. Except 
for a tiny region of the phase space, where \xik, a — u,i\ < e, such configurations can be dealt with in the 
following way: Instead of constructing pi, we construct (xik, a — Ui)Pi ■ The energy of this helper momentum 
is determined using the on-shell constraint for pi and the whole four- vector is then rescaled with l/(xifc ja — Ui). 
We typically have e of the order of 10~ 3 . 



B Enhancing photon production in the parton shower 

To improve the statistical significance of event samples with identified photons described by the parton 
shower, it is useful to enhance the corresponding branching probabilities. When doing so, one must of course 
correct for this enhancement by means of an event weight, which depends on both, acceptance and rejection 
probabilities in the parton-shower evolution. In this appendix we describe a method to incorporate an 
enhancement, which is constant over the emission phase space. Our derivation is based on the applicability 
of the veto algorithm, cf. e.g. [50] . 

Let t be the parton-shower evolution variable and f(t) the splitting kernel /C, integrated over the splitting 
variable The differential probability for generating a branching at scale t, when starting from an upper 
evolution scale t' is then given by 

V{t,t') = f(t) expj-^ df/(t)| . (17) 

A new scale t is therefore found as 

t = F- 1 [F(t')+ log R] where F(t) = J dtf(t), (18) 

and where R is a random number between zero and one. The key point of the veto algorithm is, that even if 
the integral F{t) is unknown, one can still generate events according to V using an overestimate g(t) > f(t) 
with a known integral G(t). Firstly, a value t is generated as t = G _1 [G(t') + logi?]. Secondly, the value 



2 Relation II16H is crossing invariant and can therefore be employed for all dipole configurations. 

3 For simplicity, we assume that only one splitting function exists, i.e. that there is no flavour change of the splitter during the 
evolution. The extension to flavour changing splittings is straightforward, but it would unnecessarily complicate the notation 
at this point. 
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is accepted with probaility f(t)/g(t). A splitting at t with n intermediate rejections is then produced with 
differential probability 



fit) 

V n {t,t') = J -^g{t) exp 



dig(t) I J] 



ti+1 



dti 1 



f(U) 



g(ti) exp 



dig(t) 



, (19) 



where t n+ i = t! and t§ = t. The nested integrals in Eq. (I19[) can be disentangled, and summing over n leads 
to the exponentiation of the factor g(t) — f(t), such that Eq. (fTT|) is reproduced [50] . 

Our purpose is to introduce an additional overestimate h(t) = C g(t), where C is a constant. The additional 
weight g(t)/h(t) = 1/C is then applied analytically rather than using a hit-or-miss method. This leads to 
the following expression for the differential probability to generate an emission at t with n rejections between 
t and t' 



V n (t,t') =^rh(t) exp<j- - 



d* fc© n 



1 n 



c 



gftO - /(t«)/c 



dtj 1 - 



/(*«) 



/i(ti) exp ^ - 



(20) 



The factor in the second line of Eq. (f2"U|) gives the analytic weight associated with this event, where the term 
1/C is due to the acceptance of the emission with probability f(t)/h(t). The product, which is needed for 
an exponentiation of h(t) — f(t) instead of g(t) — f(t), runs over all correction weights for rejected steps. 



C Monte-Carlo setup 

In this appendix we describe the details of the Monte-Carlo setup used to generate the results in this 
publication. We employ the SHERPA [7T] framework, which is a multi-purpose Monte-Carlo event generator 
for collider experiments. 

The matrix-element generator COMIX [81] is used to produce parton-level events for the following processes: 

e + e~ collisions fSec. HH) e+e - rr + Nr N < A max 

pp collisions (Sec. EM EM PP ~> rr + Nr N < 7V max 

The "resummed" container r implements the democratic treatment of photon and parton radiation, i.e. it 
contains the light quarks d, u, s, c and b as well as gluons and photons. In addition to all automatically 
generated tree-level amplitudes, the loop induced process gg — > 77 [3S] has been implemented. 

As a parton shower we employ the CSS [53] module which in case of N max > is merged to the matrix- 
element emissions above Q C ut = 10 GeV through the algorithm described in [62]. Unless stated otherwise, 
we employ the parton shower kinematics described in Appendix [X] All QED splitting functions are included. 
The shower cut-off has been left at the default value of p" lm = 1 GeV for the LEP runs and has been switched 
to p™ m = 2 GeV for the Tevatron runs purely for efficiency reasons. 

The PDF set employed for pp runs is CTEQ 6L (82) , which defines the corresponding a s parametrisation in 
hadron collisions. All other generator parameters are left at the default values of the Monte Carlo programs, 
since none of them is expected to have any impact on the results presented here. 

Hadron-level results for the fragmentation function analysis in Sec. 14. H are produced using the fragmentation 
module AHADIC++ [83] and the hadron and r decay package HADRONS++ [85] . The AHADIC++ default tune 
to data from the LEP experiments at the Z° resonance, obtained using the PROFESSOR [86] framework, has 
been used. To account for corrections in the ALEPH measurement, the decays of ir° and 77 have been disabled. 
Extra QED radiation in hadron decays is simulated through PHOTONS++ [ST]- All Tevatron analyses are 
presented at the parton level after parton-shower evolution. 

Multiple parton interactions have been disabled, because the presented measurements have been corrected 
for their effects. 
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